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PERGAMON 


Abstract 

A "dynamic” solar power plant (which consists of a solar collector-thermal engine 
combination) is proposed as an alternative for the more usual photovoltaic cells. A model 
for heat losses in a selective flat-plate solar collector operating on Mars is developed. An 
endoreversible Carnot cycle is used to describe heat engine operation. This provides upper 
limits for real performances. The output power is maximized. Meteorological and 
actinometric data provided by Viking Landers are used as inputs. Two strategies of 
collecting solar energy were considered: (i) horizontal collector; (ii) collector tilt and 
orientation are continuously adjusted to keep the receiving surface perpendicular on the 
Sun’s rays. The influences of climate and of various design parameters on solar collector 
heat losses, on engine output power and on the optimum sun-to-user efficiency are 
discussed. © 2000 Elsevier Science Ltd. All rights reserved. 
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1. Introduction 

A variety of future space initiatives as missions on the Moon and Mars and 
orbiting stations require the availability of power generation units with capacities 
ranging from about 1 kW to many MW utilizing chemical energy, radio-isotopes, 
solar radiation or nuclear fission as primary energy source [1,2]. The utilization of 
systems, which use internally stored nuclear and chemical energy, implies in many 
cases limited life and excessive weight. Properly designed solar power systems can 
offer long life in space. 

Among various solar technologies a leading position is played by the 
photovoltaic (PV) solar cells. However, for many years researchers realized that in 
some specific situations solar dynamic power systems with thermal energy storage 
can provide significant savings in life cycle costs and launch mass when compared 
with conventional photovoltaic power systems with battery storage [3-6]. A 
standard solar dynamic power system uses a mirror to concentrate solar radiation 
onto an absorber structure. By conduction through a solid material or circulation 
of a working fluid, the absorber heat is transferred to a thermal converter (i.e., a 
turbogenerator, Stirling engine, thermocouple or thermionic-emitter). For 
continuous operation during dark periods the use of melted materials to store 
thermal energy is being considered. Among the advantages of dynamic power 
systems one could quote their ability to provide electrical energy and heat 
simultaneously, the fact that the power plant may be unified by using either solar 
or nuclear energy or their relative invulnerability to corpuscular particles and to 
electromagnetic radiation [5]. A comparison between the various power systems 
designed for the Space Station Freedom shows the sun-to-user efficiency is 5.7- 
6.8% for the PV system and 18.7-19.8% for the dynamic system. If, in addition, 
the storage subsystem is taken into consideration, the solar dynamic system 
becomes even more attractive, mainly because the efficiency of the latent heat 
storage is higher compared to the Ni-Cd batteries of the PV system (for details 
see [7]). 

In this paper we propose a model of dynamic solar power plant operating on 
the surface of Mars. A number of power systems were already designed to be used 
under the Martian environment but all of them are based on photovoltaics. 
Indeed, PV cells are envisaged to power small rovers [8], long-endurance, remotely 
piloted aircraft capable of flight within the Martian environment [9] and to ensure 
the life support for a 40 day manned Mars surface scientific expedition [10]. More 
recently, the Sojourner of the Pathfinder mission was also powered by PV cells 
(see e.g. [11]). Our most important question is: are properly designed dynamic 
systems comparable in performance with the PV-based power systems? Some of 
the references we have already cited seem to indicate an affirmative answer. 
However, they all refer to systems operating in the interplanetary space. A 
surface-based power system has to take into consideration one very important 
feature of the Martian weather. The period between the areocentric longitude 
L s = 161-326°E and variants thereof have been referred to as the “dust storm 
season”. This period is nearly centered on perihelion, which is the time of 
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maximum insolation on Mars [12,13]. In the northern hemisphere, the Martian 
atmosphere varies greatly from year to year during the dust storm season, while at 
the other seasons most phenomena are remarkably similar from year to year. 
Zurek et al. [14] provide the most comprehensive description of Martian 
meteorology. 

The dust storm periods are characterized by an almost vanishing direct solar 
radiation (see e.g. [15] and references therein). The main consequence is that 
diffuse solar radiation cannot be concentrated and flat-plate solar collectors have 
to be used in this case. But common Earth practice proved that very low efficiency 
is associated with solar power generation using flat-plate (blackbody-like) 
collectors. The model we develop here is based on flat-plate selective solar 
collectors, which are operating at higher efficiency. 

Previous practice proved that three different cycles could be used for 
thermodynamic conversion of solar radiation: Brayton, Rankine and Stirling [3,5]. 
To provide an upper limit for the performance of an actual thermal engine we 
shall use here a Carnot endoreversible cycle. This has, on one hand, the advantage 
of generality and, on the other hand, it takes into account the irreversibilities 
associated to the heat transfer at the hot and cold parts of the engine. 


2. Solar power plant model 

Generally, a solar power plant consists of a solar collector, a thermal engine 
and a heat storage system. A solar collector mainly consists of mirror, absorber 
(receiver), transparent cover(s) and thermal insulation. The mirror is missing in 
the case of flat-plate collectors, while for some concentration systems the receiver 
is not protected by transparent cover(s). 

Direct beam (b ), diffuse (d) and ground reflected (r) fluxes of solar energy are 
incident on the collector. Their magnitude depends upon the collector’s tilt and 
orientation. Generally, the incident and absorbed solar energy density fluxes ( (p inc 


and (p a bs, respectively), are given by 


Vine = Gb + (?d + G t 

(1) 

<Pabs = P m C(rct) b G b + (m) d (G d + G r ) 

(2) 


Elere G denote solar irradiances at collector level, p m — mirror’s reflectance, C — 
concentration ratio while (ra) are effective transmittance-absorptance products. 

The solar power plant analyzed in this paper is designed to operate during the 
whole Martian year. Consequently, it consists of a flat-plate selective solar 
collector coupled to a thermal engine. No heat storage system is considered in this 
approach. One denotes by T c and T a the collector and ambient temperatures, 
respectively. 



4 


V. Badescu et al. I Renewable Energy 21 (2000) 1—22 


Solar radiation 



Fig. 1. Solar collector: (i) longitudinal section; (ii) transversal section; (iii) lateral view. 1 — absorber 
plate; 2 — carbon dioxide layer; 3 — transparent cover; 4 — bottom thermal insulation. 

2.1. Solar collector model 

A selective flat-plate solar collector is shown in Fig. 1. Its main components are: 
the absorber plate (1) with pipes for the working fluid and a transparent cover (3) 
of thickness a. The gap between (1) and (3) (thickness s) is filled with carbon 
dioxide at the Martian atmospheric pressure (carbon dioxide represents 95% of 
the Martian atmosphere). The bottom thermal insulation (4) has a thickness b. 
The collector tilt angle is (1. Various temperatures related to collector operation 
are defined in Table 1. 

The collector heat losses will be studied by using the thermal resistance method 
[2] (see Fig. 2 and Table 2). Then, the flux of heat losses from the collector 
towards the ambient is given by: 


(hot = U l (T c - 7V) = (7c - T,)/R tot 


(3) 


Table 1 

Various temperatures related to collector operation 


Temperatures 

Tc 

Absorber plate temperature 

T s i 

Temperature of lower surface of transparent cover 

T ss 

Temperature of higher surface of transparent cover 

7a 

Atmosphere temperature 

Average temperatures 

7c,si = (7c+ 7’ si )/2 

Average fluid temperature between absorber plate and lower surface of 
transparent cover (see Fig. 2) 

7’ s = (7’ si +r ss )/2 

Average transparent cover temperature 

7ss,a = (7gs 4" Ta)/2 

Average temperature of the atmosphere boundary layer near the collector 
transparent cover 

7) z = (7’ c +r a )/2 

Average bottom insulation temperature 
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Fig. 2. Solar collector thermal resistances. T a — ambient temperature; T c — absorber plate 
temperature; T s i — temperature of lower surface of transparent cover; T as — temperature of higher 
surface of transparent cover. 


where C/ L is the overall heat losses coefficient while R tot is the total thermal 
resistance. C/ L is given by (see Fig. 2): 


U L 


1 

Rrn 


1 1 
R^5 + R~6 


(4) 


where /?|_ s is computed from: 


R 1-5 = 


RiRi 
Ri + Ri 


+ -^3 + 


R 4 R 5 
R 4 + R 5 


(5) 


Details concerning the evaluation of thermal resistances follow. The thermal 
resistance R , is computed from: 


(^C0 2 }c,si^ M C,si 


Table 2 

Thermal resistances associated to solar collector operation (see Figs. 1 and 2) 


Thermal Process 

resistance 


R\ Convection and conduction between the absorbent plate (1) and the transparent 

cover (3) 

R 2 Radiation between the absorbent plate (1) and the lower surface of the transparent 

cover (3) 

R 3 Conduction through the transparent cover (3) 

R 4 Radiation between the higher surface of the transparent cover (3) and Martian 

atmosphere 

R 5 Convection between the higher surface of the transparent cover (3) and Martian 

atmosphere 

R b Conduction through the bottom thermal insulation (4) 
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Table 3 

Relationships for convection heat transfer between: (a) absorber plate and transparent cover [16]. "H' 
and ‘V' denote a horizontal and vertical solar collector; /? [rad] and L are tilt angle and collector length, 
respectively; (b) two plane parallel surfaces [17, p. 160] 


(a) 


Nuc,si = [1 — (p/2n)\Nuii + (fS/2n)Nuv 


Nu H = 0.195Gr°- 25 
Nu h = 0.068Gr°- 33 
Nuv = 0.180Gr a25 (5/L)° 11 
Nu\ = 0.065Gr a33 (j/i)° 11 


For 2 x 10 3 < Gr < 4 x 10 
For Gr>4 x 10 
For 1 x 10 4 < Gr < 2 x 10 
For 2 x 10 5 < Gr < 11 x 10‘ 


(b) 

-37/c.si — 1 

For Gr Pr < lx 10 3 

Nu c , si = 0.\05(Gr Pr) 030 

For 1 x 10 3 < Gr Pr < 1 x 10 6 

Nuc.si = 0.400(Gr Prf 20 

For 1 x 10 4 < Gr < 2 x 10 5 

Nuc M = 0.m(Gr Prf 25 

For Gr Pr > lx 10 10 


where (/co 2 )c,si is the thermal conductivity of carbon dioxide and Nu cjSi is the 
Nusselt number, both of them evaluated at the average fluid temperature 7c si . To 
compute Nuc,si, a number of criterial relationships developed for solar collectors 
were used [16] (Table 3a). For cases not covered by Table 3a, we used the general 
equations from [17] (Table 3b). 

Prandtl and Grashoff numbers in Table 3a,b were computed with: 


(fp. co,)c,sih/cxb}c,si r , guPco/( T c T si ) 

Pr = - 77 -^T---; Gr = - f — --- (7) 

WOi/C.si WC02'C,si-^C0 2 r C,si 

where g M is gravitational acceleration on Mars, while pco 2 an d Rco 2 are the 
pressure and gas constant for the carbon dioxide between absorber plate and 
transparent cover. In (7) the specific heat at constant pressure ((c P ,co 2 )c,siX the 
dynamic viscosity ((^co 2 )c,si) an d the thermal conductivity ((/hxblc.si) of carbon 
dioxide were interpolated at the average temperature 7c >si . In computations, the 
thermal properties of carbon dioxide from [8] were used. 

One neglects the radiation absorbtion in the thin carbon dioxide layer between 
absorber plate and transparent cover. The thermal resistance R 2 is then evaluated 
with: 


1/fi + I/O — 1 
2a(Tl + T%)T cm 


( 8 ) 


where and e 2 ar e the global hemispherical emittance of absorber plate and lower 
surface of transparent cover, respectively [19] and a is Boltzmann constant. 
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Table 4 

Relationships for forced convection heat transfer between transparent cover and atmosphere [17, p. 126] 


ft = 5x 10 5 ((x;co 2 >ss,a- R CO 2 7s Sja //)cO, W) 


For L < X T (laminar flow) 
For L > X T (turbulent flow) 


gO.™ p/ .0.50 

Nu ss a = 0.036 Re o so Pr 015 


The thermal resistance R 3 is given by: 


R3 = a/X s 


(9) 


The temperature dependence of transparent cover’s thermal conductivity (l s ) was 
described by using an interpolation parabolic function whose coefficients were 
obtained by using experimental data (in application the transparent cover is made 
of light silicate flint glass (65% silica wt.); thermal properties are found in [20,28]). 
The average transparent cover temperature T s was used in computations. 

The thermal resistance R 4 was evaluated with: 


R 4 = 


l/f-2 

MTfsT tWZ 


( 10 ) 


where T SSia is the average temperature of the atmosphere boundary layer near the 
collector transparent cover. 

Because of wind, the heat transfer between the transparent cover and Martian 
atmosphere is mainly by forced convection. Thus, the thermal resistance R 5 was 
computed from: 


Rs 


L 

('kxhlss.a^ss.a 


( 11 ) 


where L is collector’s length on wind direction and (Aco 2 )ss,a I s carbon dioxide 
thermal conductivity at temperature r ss>a . In (11) /V« ssai depends on the critical 
length X T which separates the laminar and turbulent wind flows over the collector. 
The relationships from Table 4 were used. In Table 4 the Reynolds and Prandtl 
numbers were evaluated with: 


_ HZ - _ (fp.COjlss.a^CxOss.a (12) 

(^/C0 2 )ss,a-^C02 ^ss,a (^C0 2 )ss,a 

where W is wind speed and the thermal properties of the atmosphere ((>?co 2 )ss,a> 
( c p, co 2 )ss,a an( l (^co 2 )ss,a) were interpolated at the average temperature T ssa . In 
applications thermal properties of carbon dioxide from [18] were used. 

The thermal resistance R 6 was evaluated with: 


Re = b/Xi- 


(13) 
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Fig. 3. Temperature — entropy diagram for an endoreversible Carnot cycle (for notations see text). 


where X iz is bottom insulation thermal conductivity, computed at the average 
temperature T iz . 

2.2. Thermal engine model 

A variety of working fluids were proposed to be used in space solar engines. 
They include mercury, potassium and rubidium in the early studies [3] or organic 
fluids and mixtures of noble gases more recently [2,4]. The latter seem to be more 
appropriate in the case of Martian solar power plants. The largest and smallest 
working fluid temperatures are denoted T' and T ", respectively (Fig. 3). 
Generally, T c > T' and T" > T. d . Flence, both 

X = T c — T' and Y = T " — T d (14) 

are positive quantities. 

The flux of solar energy ()> abs absorbed by the collector is 

</> abs = Pabs^C (15) 

where A c is the collector surface area, while cp abs is the absorbed energy density 
flux. 

Part of this flux is transferred to the working fluid ( Q c ) while the other part 
constitutes the flux of thermal losses to the ambient (g L ). Then, the solar collector 
steady-state energy balance is: 

</>abs - Q C - Qy. = 0 (16) 

The heat fluxes Qc an d Q l can be written as: 
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Q c =h c A c X and Q L = U L A C (T C - T a ) (17) 

where he and U\ are appropriately defined overall heat transfer coefficients 
between the solar collector and the working fluid and ambient, respectively. 

The Carnot engine partially converts the flux Q c into power. Thus: 

Qc - W+ Q, (18) 

where W is the output power and Q a is the thermal energy flux finally reaching 
the surroundings. There is no further increase of entropy during work production. 
Consequently 


Qc , Qa_ _ (. 

j 1 / jn 

The heat flux Q a can be written as 


(19) 


Q a — h a A a Y 


( 20 ) 


where A a and h a are appropriately defined engine’s radiator heat transfer area and 
overall heat transfer coefficient between working fluid and ambient, respectively. 


2.3. Performance indicators 


Two performance indicators were considered in this work. First, there is the 
output power W, evaluated from: 

W = h c A c X- h a A a Y (21) 


Here we used the Eqs. (17), (18) and (20). 

Second, there is the solar (or sun-to-user) efficiency defined as: 

W 

' ,solar = J~A C 

with <t>inc= (pincA C . 

Other indicators as the Carnot engine efficiency could also be used. 


( 22 ) 


2.4. Optimization 

The output power W will be maximised here. The weight of space power plants 
should be kept at minimum. Following previous studies (see e.g. [21]), we express 
the mass of the total heat transfer area (i.e. the mass M c of the solar collector and 
the mass M a of the engine’s radiator) as a k fraction from the total mass M tot of 
the space system: 


Me + M a — kM loi 


(23) 
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One denotes d c and d d the superficial mass density of collector and radiator, 
respectively. One can simply write: 

M c = d c A c ; and M d = d d A d (24) 

One can use d c = d d as a reasonable assumption. By using Eqs. (23) and (24) one 
derives the following constraint to be fulfilled: 

A c + A d = A = (25) 


Other constraints (say F\ and F 2 ) are the Eqs. (16) and (19). By using (14), (17) 
and (20) they become: 


_ hcAcX h d (A — Ac) 
1 = T c -X T d +Y 


(26) 


F i - cp dbs - U l (T c - T d ) - h c X = 0 (27) 

The power W is considered here as a function of T c , X, Y and A c . The Lagrange 
function L associated to W is: 


L = JY T ^i\F\ T ^ 2^2 


(28) 


where and X 2 a re multipliers. 

The maximum output power is obtained by solving the six equations shown in 
Table 5. Table 6 shows the solution, which applies for K > 1. There, the 
dimensionless parameters K and / are defined as: 





ffabs 

U L T d 


(29) 


Table 5 

Equations to be solved for power maximiza¬ 
tion (for details see text) 


1. 

2 . 

3. 

4. 

5. 


3 L 
3 7c 


= —3i 


hcA c X 
(T c -X f 


— X 2 UL = o 


& = hcAc + ii*^-Mc = 0 

fL=-hA A - Ac )-X 1 *gffi=0 
Mc=h C X+h lt Y+X l l 1 ^ + ^y] = 0 


dL _ hcAcX _ h a (A—Ac)Y _ /v 

dh ~~ Tc-X T a +Y ~ 


6- §£ = (Pabs ~ U L (T c ~ Ta) - h C X = 0 
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Table 6 

Optimum parameters for maximum 
output power (for details see text) 


] t _ ri t+/-(!+/ ) 1/2 ha iy 

'■ Jc.opl — 11 -I-./ Jfiy 2 (1+A - 1 / 2 ) (/ L 1 r a 

j y _ (1+/ ) — (l+/ )'^ 2 A 7 

Z - A »P' — A' '/ip+A 1/2) * c J a 


0+/ ) 1/2 -i T 
1+A 1/2 J a 


4- Tc, 0 pt = 


A 1 ' 2 
1+A '// 


A 


Table 6 and Eq. (21) allow to write the maximum power as: 

(30) 

-* a 

The optimum solar efficiency i/ so i a r,opt can be derived by using Eqs. (22), (30) and 
Table 6. 


3. Meteorologic and actinometric data 

During the application we used as input data, values of ambient pressure and 
temperature, wind speed and atmospheric optical depth measured at Viking 
Lander 1 (VLl) (22.3°N, 47.9°W) and Viking Lander 2 (VL2) (47.7°N, 225.7°W) 
sites. The data were taken from the Viking Lander Meteorology and Atmospheric 
Opacity Data Set Archive from the Planetary Data System database made 
publically available on a compact disc [22]. Three electronic files were used here, 
namely VL_MBIN.DAT and VL_TBIN.DAT (which contain values of 
atmospheric pressure and temperature and zonal and meridional wind speed, 
among others) and VL_OPAC.DAT (which contains atmospheric optical depth 
data). 

The file VL_MBIN.DAT (8325 records) contains binned and splined data 
obtained from the Viking Meteorology Instrument System (VMIS) through 
portions of the Viking Lander 1 mission. The file VL_TBIN.DAT (25750 records) 
contains binned and splined data obtained from the VMIS through most of the 
Viking Lander 2 mission and the early days of the Lander 1 mission. In both 
cases the numbering of the Martian solar days (sols) started at day 0 when each 
Lander touched down. Both files VL_MBIN.DAT and VL_TBIN.DAT contain 
mean values of meteorological parameters in each bin. The bin number starts with 
0 at local midnight and increments by 1 each 1/25 of a sol. The file 
VL_OPAC.DAT (1044 records) stores measurements of atmospheric optical depth 
and associated errors estimates. Each record contains the optical depth value and 
its associated local solar time (in Earth hours) beginning at midnight. 

As the files VL MBIN.DAT and VL TBIN.DAT are more “dense” than the 



12 


V. Badescu et al. / Renewable Energy 21 (2000) 1-22 


VL_OPAC.DAT file, a time interpolation procedure was used to obtain the values 
of the atmospheric pressure and temperature and wind velocity associated to a 
given record of atmospheric optical depth (for details see [23]). The atmospheric 
pressure values are missing for eleven records. In these cases the following 
procedure was adopted: 

1. In the case of the record VL2, year 1, areocentric longitude = 229.808 (autumn), 

sol = 200, hour =11.41, we used the arithmetic mean of the pressure values 
corresponding to sol=199, hour= 13.87 and sol = 201, hour= 11.42, 

respectively. 

2. In the case of the record VL2, year 2, areocentric longitude = 19.228 (spring), 

sol = 454, hour = 9.62, we used the arithmetic mean of the pressure values 
corresponding to sol = 453, hour= 17.33 and sol = 454, hour= 17.25, 

respectively. 

3. In the case of VL2, year 2, the first nine records during summer (areocentric 
longitudes between 93.363 and 100.510) we used the pressure value 
corresponding to the tenth record (areocentric longitude = 101.697, sol = 634, 
hour = 6.98). 

A set of complete records was finally obtained. Each complete record associates to 
a given sol number and a local solar time value a set of four meteorological 
values, namely optical depth t and atmospheric temperature T a , atmospheric 
pressure p a and wind speed W. The number of complete records is rather small 
(Table 7). 

The optical depth values were used to evaluate the incoming fluxes of solar 
radiation. The resulting file of solar energy data was used as input during the 
simulation of solar heat engine operation. The three components of global 
radiation flux (namely direct, diffuse and ground reflected) were computed by 


Table 7 

The number of complete records used during simulations of solar dynamic power system operation 




Viking Lander 1 

Viking Lander 2 

Sol numbers 

Number of 
complete records 

Sol numbers 

Number of 
complete records 

Year 1 

Spring 

- 

- 

- 

- 


Summer 

6-155 

22 

6-111 

37 


Autumn 

168-306 

108 

120-257 

77 


Winter 

307-350 

56 

269-415 

55 

Year 2 

Spring 

_ 

_ 

416-608 

147 


Summer 

- 

- 

610-784 

58 


Autumn 

- 

- 

792-872 

25 


Winter 

- 

- 

- 

- 
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using the procedure we previously proposed in [15]. The incoming fluxes of solar 
radiation were evaluated for all the complete records of Table 7. 


4. Results and discussions 

Four strategies of collecting solar energy were considered in a preliminary test 
[24]. The following two strategies will be used in this study: (i) horizontal collector 
— strategy H\ (ii) the collector tilt and orientation are continuously adjusted to 
keep the receiving surface perpendicular on Sun’s rays — strategy P. Strategy H is 
easier to use while strategy P gives the higher power output most times of the 
year. 

A number of assumptions were accepted during model application. First, a 
simple relation is used for the two effective transmittance-absorptance products: 

( ra ) b = ( Ta ) d = "tcsKps (31) 

where t cs and a ps are the transparent cover transmittance and absorber 
absorptance, respectively, both of them for short (solar) wavelengths. 

Details about solar collector design are given in Table 8. Transparent cover’s 
thermal conductivity lies between 0.691 W/rnK at — 100°C and 1.025 W/mK at 
+ 100°C [20], 

In order to allow comparison with PV power systems one can imagine a 
dynamic power system used on Mars Pathfinder’s Sojourner. Sojourner was a 
small (11.5 kg), six-wheel robotic vehicle powered by a PV array weighting 0.34 kg 
[11]. Its photovoltaic cells were designed to provide 16 W of electric power at the 
noon of a clear summer Martian day. The ratio 0.34/11.5 = 0.0296 can be used as 
a first estimate for the coefficient k to be used in Eq. (25). Preliminary tests were 
performed to determine an optimum collector area .4 C o pi close to Sojourner PV 
cell surface area (0.22 m 2 ). Finally, the value A = 0.23 m 2 was used in 
computations. It corresponds to a radiator superficial mass density d a = 1.48 kg/m 2 


Table 8 

Details about selective flat-plate solar collector design 


Quantity 

Notation 

Value 

Concentration ratio 

c 

1 

Number of transparent covers 

N 

1 

Transparent cover thickness 

a 

0.003 m 

Distance between transparent cover and absorber plate 

s 

0.045 m 

Bottom thermal insulation thickness 

b 

0.1 m 

Short wavelengths transparent cover transmittance 

T cs 

0.82 

Short wavelength transparent cover absorptance 

ftps 

0.90 

Long wavelengths (IR) absorber plate emittance 


0.10 

Long wavelengths (IR) transparent cover emittance 

£2 

0.88 

Thermal conductivity of bottom insulation (polyurethane) 


0.02 W/m K 
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Table 9 

Three design solutions for (lie system solar collector-thermal engine (/j a = heat transfer coefficient 
between the working fluid passing through the radiator of the heat engine and ambient [see Eq. (20)]; 
/?e = heat transfer coefficient between solar collector and working fluid [see Eq. (17)]) 


Case 

/; a W/m 2 K 

h c W/m 2 K 

Performance level 

a 

1 

1 

Low 

b 

10 

10 

Medium 

c 

500 

50 

High 


in Eq. (25). It is rather close to the value (4 = 2.3 kg/m 2 accepted in [25] in the 
case of the radiator of a nuclear power satellite but only half of the value d a = 3 
kg/m 2 accepted in [26] for a solar power station. One could imagine technical 
solutions decreasing the radiator thickness but keeping its mechanical rigidity (for 
instance, by using fins). 

4.1. Solar collector operation 

The solar collector thermal losses mainly depend on the heat transfer regime 
between absorber plate and transparent cover (separated through the layer 2 of 
thickness 5 in Fig. 1). We considered the case of a flat plate horizontal solar 
collector (i.e. concentration ratio C = 1). Three design solutions were analyzed 
(Table 9). They correspond to different performance levels. For analysis we have 
chosen the middle of a clear sky autumn sol, when the thermal losses are higher 
[sol 301 VF1, 14.04 local solar time (LST)]. 



Fig. 4. Overall heat loss coefficient t/ L as a function of parameter s (distance between absorber plate 
and transparent cover — see Fig. 1). Three-design solutions (a)-(c) were considered (see Table 9). 
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The dependence of the overall heat loss coefficient t/ L on the distance s between 
absorber plate and transparent cover is shown in Fig. 4 for the three cases of 
Table 9. The heat losses decrease by increasing This is specific to the conduction 
heat transfer regime. The thickness s is limited, however, by economic reasons. In 
the following we use s = 0.045 m and the high performance design solution [case 
(c) in Table 9]. Note that the chosen value of ,v is about two times larger than the 
usual value for Earth placed solar collectors. 

For Earth placed solar collectors the heat transfer between absorber plate and 
transparent cover is mainly by convection. This prompted various technical 
solutions meant to keep convection at local level (the honeycombs and Francia 
cells are examples of solutions for convection suppression (see [26] p. 404). The 
convective heat transfer is activated when the product Gr Pr exceeds 1000 (see 
Table 7). We analyzed the heat transfer above the absorber plate by using the 
whole set of available meteorological data (Table 7). 

The results show that the product Gr Pr varies between 9.65 and 211.2, well 
below the threshold value of 1000. We draw the important conclusion that on 
Mars the thermal losses between absorber plate and transparent cover are by 
conduction and radiation with convection almost negligible. 

The convection heat losses from transparent cover to atmosphere are mainly 
controlled by wind speed. The change from laminar to turbulent convection 
occurs in case the size of the heat transfer surface exceeds a certain critical length 
X T (see Table 4). The flow regime above the transparent cover was analyzed by 
using the whole set of available meteorological data. Results are shown in Fig. 5. 
The minimum value of the critical length X T is 5.43 m (Fig. 5(a)). Consequently, 
the convection heat transfer above the transparent cover is laminar. The critical 
length X T is smaller during winter and larger during summer (Fig. 5(a)). No 
significant difference between VL1 and VL2 sites is observed (note that a larger 
number of meteorological values are available for VL2 — see Table 7). The 
critical length obviously decreases when the optical depth increases (Fig. 5(b)). 
During the “clear sky” days (t < 1) the critical length values have an important 
dispersion. Scattering diminishes during the dust storm days (t > 1). 

For the design solution we selected, the values of the overall heat loss coefficient 
f/ L range between about 0.5 and 1 W/m 2 /K (Fig. 6). Broadly, these values are 
comparable with those of Earth located vacuum solar collectors (see e.g. [27], p. 
387). The dispersion of the t/ L values is slightly smaller during winter and higher 
during autumn (Fig. 6(a)). No obvious dependence of t/ L on the optical depth is 
observed (Fig. 6(b)). Flowever, the dispersion of the C/ L values is larger for values 
t < 1, i.e. during the “clear” days. Also, Fig. 6 shows that f/ L doesn’t depend 
significantly on the latitude of the solar collector. 

4.2. Thermal engine operation 

An important influence on performance is induced by the parameters h a and /? c . 
Again, a flat-plate solar collector (C = 1) with a total equivalent heat exchange 
area A = 0.23 m 2 in Eq. (25) was analyzed. We assumed s = 0.045 m. This value 
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spring: + + + VL 1 & VL2; autumn: • • • VL1; °oo VL2; 
summer: **« VL 1 & VL2; winter: □□□ VL 1; »«« VL 2; 



Areocentric longitude [ degrees ] 


(a) 


spring: * + + VL 1 & VL 2; autumn: •••VLI;ooo VL 2; 
summer: VL 1 & VL2; winter: □□□ VL 1; oo» VL2; 



Optical depth 


(b) 


Fig. 5. Critical length X T (see Table 4) as a function of (a) areocentric longitude and (b) optical depth. 
All available meteorological data were used. The design solution (c) of Table 9 was considered. 
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spring: + + + VL1 & VL 2; autumn: ••• VL1; °°o VL 2; 
summer: »« VL 1 & VL 2; winter: »»" VL 1 ; ooo VL 2; 



(a) 


spring: ** + VL 1 & VL2; autumn: • • • VL 1; ooo VL2; 
summer: VL 1 & VL2; winter :VL 1 ; oo® VL2; 



(b) 


Fig. 6. Overall heat loss coefficient [/ L as a function of (a) areocentric longitude and (b) optical depth. 
All available meteorological data were used. The design solution (c) of Table 9 was considered. 
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K 

(a) 



(b) 

Fig. 7. Maximum power W max provided by the thermal engine as a function of the overall heat transfer 
coefficients /j a and h c . (a) horizontal solar collector (strategy H); (b) solar collector permanently 
perpendicular on the Sun's rays (strategy P). Unconcentrated solar radiation was considered (VLl, 
Autumn, sol 301, 14.04 LST). 
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h a 

(a) 



h a 

(b) 


Fig. 8. Solar efficiency (? so lar.opt as a function of the overall heat transfer coefficients /; a and he- (a) 
horizontal solar collector (strategy H)\ (b) solar collector permanently perpendicular on the Sun’s rays 
(strategy P ). Unconcentrated solar radiation was considered (VL1, Autumn, sol 301, 14.04 LST). 
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is a compromise between a good thermal insulation and a reasonable size (see 
Fig. 4). For analysis we have chosen the middle of a clear sky autumn sol [sol 301 
VL1, 14.04 local solar time (LST)]. Some results are shown in Figs. 7 and 8 for 
the maximum power provided by the thermal engine and the solar efficiency, 
respectively. 

Generally, the influence on performance is strong for low values of both heat 
transfer coefficients h a and he (less than 100 W/m 2 /K). The performance is 
practically not dependent on the higher values of h c . The influence of h a on 
performance is more significant. 

The strategies of collecting solar radiation have a rather weak influence on the 
dependence of performance on h a and h c (see Fig. 7(a) and (b), on the one hand, 
and Fig. 8(a) and (b), on the other hand). As expected, strategy P leads to higher 
performance. 

The maximum output power could be as high as 12 W and 13 W in case of 
strategies H and P of collecting solar radiation, respectively. This is lower than 
the designed power of Sojourner’s PV array (16 W at noon). FLowever, this 
comparison put in disadvantage the solar dynamic system as we used input 
meteorological data from late autumn at VL1 site while Sojourner was designed to 
operate during Martian summer. This is a consequence of the fact that none of 
the 22 complete records available during summer at VL1 site (Table 7) 
corresponds to solar noon. 

The influence of both h a and h c is more important on the maximum power 
provided by the heat engine than on the overall efficiency (compare, for instance, 
Figs. 7(a) and 8(a)). 


5. Conclusions 

The solar power plant analyzed in this paper is designed to operate during the 
whole Martian year. Consequently, it consists of a selective flat-plate solar 
collector coupled to a thermal engine operating on an endoreversible Carnot cycle. 
Methods of finite time thermodynamics are used to maximize the output power. 
During the optimization procedure the weight of the space power plant was used 
as a constraint. 

During the application we used as input data values of atmospheric pressure 
and temperature, wind speed and atmospheric optical depth measured at Viking 
Landers sites. 

The following two strategies of collecting solar energy were analyzed in this 
study: (i) horizontal collector — strategy H\ (ii) the collector tilt and orientation 
are continuously adjusted to keep the receiving surface perpendicular on Sun’s 
rays — strategy P. 

The main results reported in this work are as follows. 

1. The heat losses between the solar collector’s absorber plate and transparent 
cover are by conduction and radiation. This differs from the Earth-based solar 
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collector where the main heat loss mechanism is convection. 

2. The convection heat transfer above the transparent cover is laminar. This 
differs from the Earth-based solar collectors where turbulence is the main 
convection mechanism. 

3. For the design solution we selected, the values of the overall heat loss 
coefficient t/ L range between about 0.5 and 1 W/m 2 /K. Roughly, these values 
are comparable with those of Earth located vacuum solar collectors. The 
dispersion of the C/ L values is slightly smaller during winter and higher during 
autumn. No obvious dependence of t/ L on the optical depth is observed. 
Elowever, the dispersion of the C/ L values is larger for values during the “clear” 
days. t/ L doesn’t depend significantly on the latitude of the solar collector. 

4. Two parameters are used to quantify the heat transfer inside the thermal 
engine. They are h a ( = heat transfer coefficient between the working fluid 
passing through the radiator of the heat engine and ambient; see Eq. (20)) and 
h c (= heat transfer coefficient between solar collector and working fluid [see Eq. 
(17)]). The influence of h a on performance is more significant. 

5. The strategy of collecting solar radiation has a rather weak influence on the 
dependence of performances on h a and h c . As expected, strategy P leads to 
higher performances. 

6. The dynamic solar power plants equipped with selective flat-plate collectors 
could provide power comparable to that of PV cell systems. 
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